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Parametric  Adjustments 
to  the  Rankine  Vortex  Wind 
Model  for  Gulf  of  Mexico 
Hurricanes 

Parametric  wind  models  are  often  used  to  reconstruct  hurricane  wind  fields  from  a  lim¬ 
ited  set  of  hurricane  parameters.  Application  of  the  Rankine  Vortex  and  other  models 
used  in  forecasting  Gulf  of  Mexico  hurricanes  show  considerable  differences  between  the 
resulting  wind  speeds  and  data.  The  differences  are  used  to  guide  the  development  of 
adjustment  factors  to  improve  the  wind  fields  resulting  from  the  Rankine  Vortex  model. 
The  corrected  model  shows  a  significant  improvement  in  the  shape,  size,  and  wind  speed 
contours  for  14  out  of  17  hurricanes  examined.  The  effect  on  wave  fields  resulting  from 
the  original  and  modified  wind  fields  are  on  the  order  of  4  m,  which  is  important  for  the 
estimation  of  extreme  wave  statistics.  [DOI:  10.1115/1.4006148] 


1  Introduction 

In  the  Gulf  of  Mexico  (GOM),  extreme  wave  height  estimates 
used  until  recently  for  the  design  of  offshore  structures  were  pro¬ 
vided  by  the  API  [1,2].  The  recommended  100-yr  design  signifi¬ 
cant  wave  height  (SWH)  has  been,  until  1990,  of  the  order  of  11 
m.  Comparable  estimates  of  the  100-yr  conditions  were  obtained 
by  Palao  et  al.  [3]  and  Panchang  et  al.  [4],  who  have  reported  con¬ 
siderable  spatial  variability  in  these  estimates.  However,  these 
estimates  have  been  exceeded  since  2004  by  several  meters.  Nota¬ 
ble  examples  include  Hurricanes  Ivan,  Katrina,  and  Rita  where 
SWHs  exceeding  16  m  were  recorded  at  various  locations  (e.g., 
Jeong  and  Panchang  [5];  Panchang  and  Li  [6]).  At  other  locations, 
smaller  SWHs  were  recorded,  but  they  still  exceeded  the  100-year 
estimates  by  several  meters.  As  a  result,  it  has  become  necessary 
to  reassess  the  extreme  wave  climatology  in  the  GOM  region 
(e.g.,  Berek  et  al.  [7]).  Similar  efforts  are  underway  to  reassess 
extreme  storm  surge  and  coastal  inundation  effects  in  this  region. 

The  U.S.  Army  Corps  of  Engineers,  Geological  Survey,  and 
National  Oceanic  and  Atmospheric  Administration  (NOAA)  are 
developing  methods  to  predict  the  coastal  impacts  of  extreme 
storms  on  the  coasts  of  the  United  States.  The  primary  emphasis 
in  these  studies  is  on  the  numerical  modeling  of  hurricanes 
(https://ipet.wes.army.mil/)  and  typhoons  [8-12],  whereas  Demi¬ 
rbilek  et  al.  [13]  and  Demirbilek  and  Nwogu  [14]  focused  on  mod¬ 
eling  the  effects  of  highly  nonlinear  storm  waves  on  fringing 
reefs.  These  efforts  are  directed  towards  predicting  where  maxi¬ 
mum  coastal  erosion  will  occur,  where  storm  surge  and  waves 
will  overtop  beaches,  where  sand  dunes  will  retreat  landward,  and 
where  breaches  will  sever  barrier  islands  and  create  new  inlets. 
For  example,  during  Hurricane  Katrina,  the  surge  level  in  Lake 
Pontchartrain  was  roughly  the  same  as  the  design  levels  assumed 
for  the  Hurricane  Protection  System.  On  the  east  side  of  New  Or¬ 
leans,  Katrina-generated  surges  were  significantly  greater  than  the 
design  criteria,  ranging  from  5.2  to  6.1  m  compared  to  the  3.7  to 
4.3  m  assumed  in  the  design.  Typically,  these  types  of  studies 
require  the  development  of  historical  wind  fields,  which  can  then 
be  used  for  developing  design  criteria  for  offshore  structures  and 
coastal  protection  systems  such  as  levees. 

In  the  GOM,  wind/wave  data  can  be  obtained  from  NDBC 
buoys  which  provide  data  for  a  maximum  duration  of  approxi- 
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mately  35  years.  These  data  provide  “spot  measurements”  and  at 
other  locations  one  must  resort  to  numerical  modeling  to  obtain 
the  correct  spatial  variability  in  the  wind  (and  in  the  resulting 
wave  and  storm  surge)  estimates.  For  numerical  modeling  pur¬ 
poses,  four  types  of  wind  fields  (on  different  spatial  and  temporal 
grids)  are  available: 

(1)  Using  the  most  sophisticated  models  available,  the  National 
Center  for  Environmental  Prediction  (NCEP)  produces 
wind  fields  on  an  ongoing  basis  every  6  h.  In  the  GOM,  the 
simulations  are  made  using  their  “Western  North  Atlantic” 
and  “North  Atlantic  Hurricane”  models,  and  wind  fields  are 
available  on  a  0.25°  x  0.25°  grids.  These  simulations  repre¬ 
sent  the  estimates  of  the  entire  GOM  wind  fields;  however, 
they  are  not  available  prior  to  1999.  Thus,  alternative  meth¬ 
ods  must  be  explored  if  hindcasts  for  that  period  are  needed. 

(2)  The  NCEP  and  the  National  Center  for  Atmospheric 
Research  have  developed  the  “Reanalysis”  wind  field  data¬ 
set,  using  a  combination  of  mathematical  models  and  data 
assimilation  (Kalnay  et  al.  [15]).  An  example  is  shown  in 
Fig.  1  (left)  for  Hurricane  Gordon.  These  wind  fields  are 
available  for  the  period  starting  at  1948  at  a  temporal  reso¬ 
lution  of  6  h;  however  the  spatial  resolution  is  coarse 
(2.5°  x  2.5°).  As  a  result,  some  features  of  a  hurricane  may 
not  be  well  represented  by  these  data,  despite  their  use  in 
large  areas  such  as  the  Atlantic  [16-20]. 

(3)  A  dataset  representing  hurricane  measurements  since  1995 
has  been  developed  by  the  National  Hurricane  Center  [21,22]. 
This  dataset,  called  H*Wind,  is  available  for  the  post- 1994  pe¬ 
riod,  and  has  been  widely  used  by  researchers  for  various 
applications  (e.g.,  Kennedy  et  al.  [23];  Powell  et  al.  [24]).  It  is 
an  estimate  of  the  wind  field  based  on  available  observations, 
viz.,  aircraft-based,  land-based,  sea-based,  and  satellite-based. 
Based  on  a  standardization  technique  to  process  data  from 
diverse  sources,  it  provides  wind  fields  at  a  resolution  of 
approximately  6  km  (an  example  is  shown  in  Fig.  1).  As  may 
be  expected,  however,  this  dataset  is  not  continuous  (except  in 
the  recent  past),  does  not  cover  all  hurricanes,  exhibits  a  spa¬ 
tial  range  limited  to  the  immediate  vicinity  of  the  hurricane, 
and  is  available  at  irregular  time  steps. 

(4)  A  dataset  pertaining  to  a  limited  set  of  storm  parameters 
has  been  developed  by  the  NOAA  for  a  period  going  back 
to  the  second  half  of  the  19th  century.  This  dataset,  called 
HURDAT,  provides,  at  6  h  intervals,  information  such  as 
the  location  of  the  storm  center  (LatC,  LonC),  storm  direc¬ 
tion  {Qs),  storm  speed  (U5),  maximum  wind  speed  ( Vm ),  and 


NOVEMBER  2012,  Vol.  134  /  041102-1 

2012  by  ASME 


Downloaded  31  Jul  2012  to  140.194.40.43.  Redistribution  subject  to  ASME  license  or  copyright;  see  http://www.asme.org/terms/Terms_Use.cfm 


Report  Documentation  Page 

Form  Approved 

OMB  No.  0704-0188 

Public  reporting  burden  for  the  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources,  gathering  and 
maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of  information, 
including  suggestions  for  reducing  this  burden,  to  Washington  Headquarters  Services,  Directorate  for  Information  Operations  and  Reports,  1215  Jefferson  Davis  Highway,  Suite  1204,  Arlington 

VA  22202-4302.  Respondents  should  be  aware  that  notwithstanding  any  other  provision  of  law,  no  person  shall  be  subject  to  a  penalty  for  failing  to  comply  with  a  collection  of  information  if  it 
does  not  display  a  currently  valid  OMB  control  number. 

1.  REPORT  DATE 

NOV  2012  1 REPORT  TYPE 

3.  DATES  COVERED 

00-00-2012  to  00-00-2012 

4.  TITLE  AND  SUBTITLE 

Parametric  Adjustments  to  the  Rankine  Vortex  Wind  Model  for  Gulf  of 
Mexico  Hurricanes 

5a.  CONTRACT  NUMBER 

5b.  GRANT  NUMBER 

5c.  PROGRAM  ELEMENT  NUMBER 

6.  AUTHOR(S) 

5d.  PROJECT  NUMBER 

5e.  TASK  NUMBER 

5f.  WORK  UNIT  NUMBER 

7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

U.S.  Army  Engineer  Research  and  Development  Center, Coastal  and 
Hydraulics  Laboratory ,3909  Halls  Ferry  Road, Vicksburg, MS, 39180 

8.  PERFORMING  ORGANIZATION 

REPORT  NUMBER 

9.  SPONSORING/MONITORING  AGENCY  NAME(S)  AND  ADDRESS (ES) 

10.  SPONSOR/MONITOR’S  ACRONYM(S) 

11.  SPONSOR/MONITOR’S  REPORT 
NUMBER(S) 

12.  DISTRIBUTION/AVAILABILITY  STATEMENT 

Approved  for  public  release;  distribution  unlimited 

13.  SUPPLEMENTARY  NOTES 

14.  ABSTRACT 

15.  SUBJECT  TERMS 

16.  SECURITY  CLASSIFICATION  OF:  17.  LIMITATION  OF 

ARSTRATT 

1 8 .  NUMBER  1 9a.  NAME  OF 

OF  PAGES  RESPONSIBLE  PERSON 

a.  REPORT  b.  ABSTRACT  c.  THIS  PAGE  Same  aS 

unclassified  unclassified  unclassified  Report  (SAR) 

17 

Standard  Form  298  (Rev.  8-98) 

Prescribed  by  ANSI  Std  Z39-18 


Longitude  (W) 


Fig.  1  Wind-field  (m/s)  for  Hurricane  Gordon,  0000  UTC  17  Sept.  2000;  Reanalysis  (left),  and  H*Wind  field  (right) 


storm  central  pressure  (Pc).  Obviously,  these  data  do  not 
provide  details  of  a  complete  wind  field. 

For  recreating  historical  wind  fields  on  a  continuous  basis  for  a 
grid  covering  the  entire  GOM,  neither  the  H*Wind  nor  the  Rean¬ 
alysis  wind  fields  are  completely  adequate.  Planetary  boundary 
layer  models  may  be  used  to  construct  spatially  varying  wind 
fields  for  the  GOM  region  (e.g.,  Thompson  and  Cardone  [12]),  but 
this  requires  using  a  complicated  wind  model  that  may  not  be 
practical.  Therefore,  here  we  have  taken  the  approach  of  develop¬ 
ing  wind  fields  using  parametric  formulations  that  depend  on  a 
limited  set  of  parameters  (HURD AT).  This  approach  has  been 
used,  for  instance,  by  Phadke  et  al.  [25]  and  Sanchez  et  al.  [9,10], 
for  simulating  tropical  cyclones  near  Hawaii,  and  by  McAfee  and 
Pearson  [26]  for  simulating  mid- Atlantic  hurricanes.  Additionally, 
extensions  to  real-time  surface  wind  forecasting  using  such  mod¬ 
els  is  possible  (Xie  at  al.  [27]). 

A  summary  of  five  parametric  models  is  given  by  MacAfee  and 
Pearson  [26].  These  are  the  Rankine  Vortex  (RV)  model  [25],  the 
SLOSH  model  [28],  the  Holland  model  [29],  the  vortex  simulation 
model  [30],  and  the  Willoughby  and  Rahn  model  [31].  Sanchez 
et  al.  [9,10]  have  used  variations  of  these  parameteric  models  in 
their  studies  of  typhoon-induced  flooding  and  inundation  of  the  Pa¬ 
cific  Islands  by  tropical  storms,  and  provided  corrections  for  asym¬ 
metric  wind  patterns  and  options  for  the  parameterization  of  storm 
characteristics.  In  addition,  Emanuel  et  al.  [32]  describe  an  alterna¬ 
tive  model  with  a  structural  form  similar  to  the  parametric  vortex 
developed  by  Holland  [29].  In  general,  these  parametric  models 
provide  the  wind  speed  V(x,  y )  as  a  function  of  hurricane  parame¬ 
ters  Vm  and  Pc,  maximum  hurricane  wind  speed  and  central  pres¬ 
sure,  and  also  other  model-calculated  quantities  such  as  the  radius 
to  maximum  wind  (Rm),  sea  level  pressure  at  the  last  closed  isobar 
( Pn ),  etc.  (Note  that  all  numerical  models  do  not  include  all  these 
parameters.)  As  an  example,  for  Hurricane  Katrina,  HURDAT  pro¬ 
vides  the  parameters  noted  at  the  bottom  of  this  page. 

These  parameters,  in  combination  with  various  parametric  mod¬ 
els,  result  in  the  wind  field  plots  shown  in  Fig.  2  (left).  Typically, 
the  basic  models  yield  a  symmetric  pattern  and  asymmetry  may 


also  be  parametrically  introduced  [33],  which  leads  to  the  plots 
shown  in  Fig.  2  (right). 

There  is,  of  course,  no  consensus  on  which  model  should  be 
used;  nor  is  it  reasonable  to  expect  any  one  model  to  faithfully 
reproduce  every  hurricane.  Researchers  have  had  to  make  adjust¬ 
ments  to  suit  individual  cases  or  locations,  as  is  noted  later.  For 
the  purpose  of  recreating  historical  wind  fields  in  the  GOM,  here 
we  focus  on  the  RV  model  (described  in  Sec.  2).  This  choice  is 
based  on  the  conclusions  of  Phadke  et  al.  [25]  who  found  that  the 
first  three  models  produced  the  same  results  in  the  core  of  the  hur¬ 
ricane,  however,  the  RV  model  had  better  overall  agreement  with 
data  outside  the  core.  The  performance  of  this  model  is  examined 
against  data  pertaining  to  five  recent  hurricanes  in  Sec.  2.  The 
observed  discrepancies  are  used  to  make  adjustments  to  the  RV 
model;  these  are  described  in  Sec.  3.  In  Sec.  4,  the  performance  of 
the  adjusted  model  is  examined  using  wind  data  from  twelve  other 
hurricanes.  The  effect  of  the  modified  wind  fields  on  wave  simula¬ 
tion  is  also  demonstrated. 


2  The  RV  Model  and  its  Application  to  Recent 
Hurricanes 


The  RV  model  provides  a  radially  symmetric  hurricane  wind 
field  as  follows 


V  = 


R  <  Rm 
R  >  Rm 


(1) 


where  Vm  is  the  maximum  wind  speed,  R  is  the  radial  distance 
from  the  storm  center,  Rm  is  the  radius  to  maximum  winds,  and  B 
is  a  shape  parameter  (?^0.5).  The  radius  to  maximum  winds  ( Rm 
in  km)  has  been  described  by  MacAfee  and  Pearson  [26]  as  a 
function  of  the  latitude  (cp  in  degrees),  the  central  pressure  (mb), 
and  the  pressure  along  the  last  closed  isobar  (Pn  in  mb)  as  follows 

In Rm  =  2.636  -  0.00005086(P„  -  Pc)2+  0.0394899 <p  (2) 


Translational 

Max  wind  speed: 

Pressure: 

Month 

Day 

Hour 

LatC 

LogC 

Dir.  ((h) 

speed:  Vs  (km/h) 

Vm  (km/h) 

Pc  (mb) 

Type 

8 

28 

12 

25.7 

87.7 

300 

18 

270 

909 

Category  5 
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Fig.  2  Symmetric  (left)  and  asymmetric  (right)  wind  speeds  (m/s)  during  Hurricane  Katrina  (1200  UTC  28  Aug. 
2005);  (a)  Rankine  Vortex  model,  ( b )  SLOSH  model,  and  (c)  Holland  model 


where 


Pn  =  Pc  -  20.69  +  1.33Vm  +  0.11<p 


(3) 


Vrv=  {  (-Lcos  0r)2+ 


V  sin  6r  + 


VsRmR 


1  2  ' 


(R2m+R2) 


1/2 


(4) 


Willoughby  et  al.  [34]  provide  an  alternative  formula  to  estimate 
Rm  as  a  function  of  the  latitude  and  maximum  wind  speed. 
The  importance  of  properly  introducing  asymmetry  to  the  RV 
model  has  been  emphasized  by  Liu  et  al.  [35].  This  can  be 
accomplished  by  making  an  adjustment  to  V  due  to  storm  trans¬ 
lation.  The  net  motion-adjusted  wind  velocity  may  then  be 
obtained  as 


where  9r  is  the  angle  between  the  storm  direction  and  the  radius 
to  a  particular  grid  point,  and  V  is  the  estimated  model  wind 
speed,  resulting  from  Eq.  (1)  [26].  For  the  wind  direction,  a  circu¬ 
lar  wind  flow  pattern  is  assumed  and  a  cross-isobaric  direction 
correction  /?,  proposed  by  Bretschneider  [36],  can  be  applied  to 
each  grid  point 
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Table  1  HURDAT  dataset  for  five  hurricanes 
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12 
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86 
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18 
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14 
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9 
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9 

15 

6 
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16 
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9 

15 

12 
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340 

20 
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9 

15 

18 
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22 
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23 

18 

27.8 

92.3 

315 

16 

205 

930 

Category  3 

9 

24 

0 

28.6 

93 

320 

18 

195 

931 

Category  3 

DOLLY  -  2008 

Translational 

Max  wind 

Pressure 

Month 

Day 

Hour 

Lat. 

Log. 

Dir. 

speed  (km/h) 

speed  (km/h) 

(mb) 

Type 

7 

21 

18 

22.8 

90.4 

305 

31 

85 

1005 

Tropical  Storm 

7 

22 

0 

23 

92 

280 

25 

85 

1000 

Tropical  Storm 

7 

22 

6 

23.2 

93.3 

280 

22 

85 

999 

Tropical  Storm 

7 

22 

12 

23.7 

94.1 

305 

14 

100 

993 

Tropical  Storm 

7 

22 

18 

24.3 

94.9 

310 

16 

110 

990 

Tropical  Storm 

7 

23 

0 

24.9 

95.7 

310 

16 

120 

982 

Category  1 

7 

23 

6 

25.4 

96.2 

320 

11 

130 

982 

Category  1 

IKE  -  2008 

Translational 

Max  wind 

Pressure 

Month 

Day 

Hour 

Lat. 

Log. 

Dir. 

speed  (km/h) 

speed  (km/h) 

(mb) 

Type 

9 

10 

0 

23.1 

84 

300 

12 

120 

968 

Category  1 

9 

10 

6 

23.4 

84.6 

300 

11 

130 

964 

Category  1 

9 

10 

12 

23.8 

85.2 

305 

11 

150 

959 

Category  1 
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Table  1 .  Continued 


IKE  -  2008 


Month 

Day 

Hour 

Lat. 

Log. 

Dir. 

Translational 
speed  (km/h) 

Max  wind 
speed  (km/h) 

Pressure 

(mb) 

Type 

9 

10 

18 

24.2 

85.8 

305 

11 

160 

958 

Category  2 

9 

11 

0 

24.7 

86.4 

310 

12 

160 

944 

Category  2 

9 

11 

6 

25.1 

87.1 

300 

12 

160 

945 

Category  2 

9 

11 

12 

25.5 

88 

295 

16 

160 

946 

Category  2 

9 

11 

18 

25.8 

88.9 

290 

14 

160 

952 

Category  2 

9 

12 

0 

26.1 

90 

285 

18 

160 

954 

Category  2 

9 

12 

6 

26.4 

91.1 

285 

18 

165 

954 

Category  2 

9 

12 

12 

26.9 

92.2 

295 

20 

175 

954 

Category  2 

9 

12 

18 

27.5 

93.2 

305 

18 

175 

954 

Category  2 

9 

13 

0 

28.3 

94 

320 

18 

175 

952 

Category  2 

{10°  +  ^1  +  ^ ,  R  <  Rm 

( R  \  (5) 

20°  +  25°  f  —  -  1  j ,  Rm<R<  1.2 Rm 

25°,  m  R>  1.2 Rm 

The  final  wind  flow  is,  thus,  anticlockwise  at  a  direction  of  (90-/3) 
degrees  to  the  radius  at  any  point. 

In  this  study,  five  prominent  Gulf  of  Mexico  hurricanes  that 
occurred  during  the  2004-2008  period  are  chosen  for  a  detailed 
examination.  These  are  Ivan  (2004),  Katrina  (2005),  Rita  (2005), 
Dolly  (2008),  and  Ike  (2008).  At  the  location  of  one  or  more 
NDBC  buoys  in  the  GOM,  wave  heights  recorded  during  four  of 
these  storms  were  larger  than  those  recorded  prior  to  2004,  and 
in  some  cases  exceeded  the  100-yr  event  [5,6].  The  HURD  AT 
details  of  these  hurricanes  are  provided  in  Table  1.  As  can  be 
seen,  these  data  provide  a  total  of  56  snapshots.  (In  actuality, 
more  snapshots  are  available  in  the  HURDAT  dataset,  but  the 
others  extend  over  regions  outside  the  GOM,  including  land,  and 
hence,  are  not  considered.)  These  snapshots  cover  the  entire 
range  of  hurricane  development  phases/levels  (e.g.,  from  the 
“tropical  storm”  category  to  category  5,  central  pressures  rang¬ 
ing  between  902  mb  and  1005  mb,  and  wind  speeds  varying 
between  85  km/h  and  285  km/h).  Additionally,  these  storms 
were  chosen  because  they  had  the  largest  number  of  HURDAT 
snapshots  for  validation. 

Some  results  of  the  wind  fields  calculated  using  Eqs.  (l)-(5)  for 
the  set  of  parameters  in  Table  1  are  shown  Figs.  3(a)  and  4(a). 
The  results  represent  one  snapshot  from  each  of  the  five  storms.  The 
corresponding  H*Wind  datasets  are  shown  in  Figs.  3(b)  and  4(b).  A 
comparison  suggests  that  the  parametric  model  needs  adjustment  to 
obtain  a  better  match  with  the  data.  In  particular,  for  some  storms  the 
“bean  shape”  in  the  core  is  insufficiently  developed  and  the  distribu¬ 
tion  of  the  velocity  contours  is  tighter  (smaller)  than  the  data  (e.g., 
the  plots  show  considerable  underestimation  in  many  areas,  espe¬ 
cially  on  the  right  side  of  the  hurricane).  Figure  5  shows  the  wind 
speed  along  the  central  East-West  radial,  calculated  using  not  only 
the  RV  model  (Eqs.  (l)-(5)),  but  also  the  widely-used  SLOSH  and 
Holland  models.  All  three  models  show  a  mismatch  (and  a  similar 
mismatch  was  also  seen  in  the  N-S  transect)  indicating  a  need  for 
adjustment.  Wind  speed  errors  are  as  large  as  15  m/s,  and  when  used 
with  a  wave  model,  the  resulting  errors  in  SWHs  can  be  of  the  order 
of  several  meters,  as  shown  in  Fig.  6. 

This  type  of  mismatch  has  also  been  observed  in  other  studies 
and  typically,  adjustments  may  have  to  be  made  by  using  data  to 
develop  alternative  equations.  For  example,  Xie  et  al.  [27]  made 
adjustments  to  the  Holland  model  based  on  four  mid- Atlantic  hur¬ 
ricanes.  Instead  of  treating  Rm  as  independent  of  6,  they  found  it 
preferable  to  describe  it  as  a  power  series  of  6  with  the  constants 
being  determined  on  a  case-by-case  basis,  for  which  details  such 
as  wind  speeds  at  specific  locations  are  required.  For  our  work, 
such  data  are  not  available  (especially  for  the  pre-1994  period). 
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Sanchez  et  al.  [9,10],  MacAfee  and  Pearson  [26],  and  Willoughby 
et  al.  [34]  also  made  adjustments  which  were  tailored  for  mid¬ 
latitude  applications. 

3  Adjustments  to  the  RV  Model 

In  order  to  make  adjustments  to  the  RV  model,  the  calculated 
wind  speed  plots  were  compared  with  H*Wind  plots  for  the  five 
chosen  hurricanes.  In  general,  the  RV-model  results  were  fre¬ 
quently  similar  to  H*Wind  for  strong  hurricanes  (e.g.,  for 
Pc  <  930  mb).  In  other  instances,  however,  both  the  shape  and  the 
magnitudes  were  incorrect.  Based  on  the  differences,  correction 
factors  d(6)  and  Cb  were  established.  The  first  correction  factor 
was  intended  to  modify  the  RV-modeled  hurricane  shape  and  the 
second  was  intended  to  modify  the  size  and  velocity  distributions. 
The  net  wind  field  VCrv  is  obtained  in  two  stages,  as  follows 

VRVx  =  Vgy  x  d(6) 

Vcrv  =  Vrv i  x  Cb 

3.1  Development  of  Correction  Factor  d(0).  A  detailed  ex¬ 
amination  of  the  modeled  and  H*Wind  plots  indicated  that  the 
comparison  was  acceptable  for  lower  Pc  values  (Fig.  7  (top)), 
however,  for  larger  values  of  Pc  the  hurricane  shape  was  more 
deformed  relative  to  the  modeled  shape  (Fig.  7  (bottom)),  i.e.,  the 
hurricane  was  less  circular  in  shape  and  developed  a  more  pro¬ 
nounced  bean-shaped  central  feature  than  in  the  model. 

To  accommodate  these  features,  it  was  necessary  to  appropri¬ 
ately  deform  the  modeled  RV  contour  shapes.  Since  the  RV  model 
produces  largely  circular  contours,  the  values  of  the  H*Wind  data 
velocities  along  any  circle  (after  appropriate  normalization)  can 
be  used  as  a  measure  of  the  required  deformation  d(6).  For  exam¬ 
ple,  denoting  the  H*Wind  velocities  at  three  arbitrary  points  P,  Q, 
and  B  along  a  circle  (shown  in  Fig.  8)  by  VP ,  Vq ,  and  VB,  respec¬ 
tively,  the  ratios  dP  —  VP/VB  and  dQ  =  V q/V B  can  be  used  as  a 
measure  of  the  deviation  of  the  contour  from  a  circle.  A  circle 
with  radius  4 Rm  is  nominally  selected  for  this  purpose  (to  go  suffi¬ 
ciently  far  from  the  center),  and  d(0)  values  are  estimated  at 
selected  points  on  the  circle.  These  estimates  are  used  to  calculate 
the  deformation  d(6)  at  other  points  by  curve-fitting. 

To  be  specific,  the  point  with  the  smallest  velocity  on  this  circle 
was  denoted  by  “C.”  Two  other  points  (A  and  B)  were  located  on 
the  circle  in  order  to  create  radii  perpendicular  to  the  one  to  C 
(Fig.  8  (left)).  Using  point  B  as  the  base  velocity,  the  ratios 
dA  =  VA/VB  and  dc  =  VC/VB  were  estimated  to  provide  measures 
of  deformation  at  these  points.  These  data  are  provided  in  Table  2 
in  the  form  of  averages  for  each  storm. 

The  values  of  dA  for  the  56  plots  were  found  to  vary  in  a  rela¬ 
tively  small  range  (±4%)  around  the  average,  and  hence,  the  aver¬ 
age  (?^1.15)  was  used  as  the  deformation  factor.  On  the  contrary, 
the  ratio  dc  displayed  less  uniformity  (varying  between  —16% 
and  +6%  around  the  average)  and  was  found  to  depend  on  the 
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(a) 


(b) 


Fig.  3  Wind  fields  (m/s)  for  hurricanes  Ivan,  Katrina,  and  Rita:  (a)  asymmetric  RV  model,  and  ( b )  H*Wind 


central  pressure,  as  seen  in  Fig.  9.  Since  930  <  Pc  <  970  for  most 
of  the  cases  examined,  a  best-fit  line  was  used  to  define  this  defor¬ 
mation  ratio  as  follows 

dc  =  0.8  +  02. a  for  930  <PC<  970 

dc  =  1,  for  Pc  <  930  (6) 

dc  =  0.8,  for  Pc  >  970 

041102-6  /  Vol.  134,  NOVEMBER  2012 


where  a  =  (970  —  Pc)  /  40. 

It  is  now  necessary  to  find  d{6)  using  the  deformation  at  A,  B, 
and  C(dA ,  dB,  and  dc),  for  which  reference  to  Fig.  10  is  made. 
The  right  curve  between  A  and  B  was  described  by 
d(0)  =  ai0  +  bi,  where  dA  and  dB  (i.e.,  information  at  6  =  0  and 
6  =  n)  are  used  to  determine  ax  and  bx.  The  left  half  was  described 
by  d(6)  =  a2  cos(2 6 )  +  b2  for  n<q<  3/2 p  and  d(6)  =  a3 
cos(2 0)  +  b3  for  3/271  <q<2p,  and  a2  and  b2  are  obtained  from 
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(a)  (b) 


Fig.  4  Wind  fields  (m/s)  for  hurricanes  Dolly  and  Ike:  (a)  asymmetric  RV  model,  and  ( b )  H*Wind 


dB  and  dc,  and  a3,  b3  are  from  dc  and  dA.  These  calculations  lead 
to 

d(0)  =  1.15  -  (O.15/7r)0,  0  <  0  <  n 

d(6)  =  1  —  (1  —  dc)  sin2  0,  71  <  0  <  ^71  ^ 

3 

d(Q)  =  dc  +  (1.15  —  dc)  cos2  0,  -n  <  0  <2n 

Multiplying  the  RV  model  velocities  by  d(6)  changes  the  largely 
circular  contour  shape;  however,  it  is  first  necessary  to  determine 
the  locations  with  the  closest  correspondence  to  A,  B,  and  C  on 
the  modeled  contours.  For  the  parametric  model,  theoretically 
(based  on  Eq.  (4))  the  radius  to  point  A  is  in  the  same  direction  as 
the  storm  direction.  However,  the  data  show  that  this  is  not  always 
true.  The  actual  storm  direction  can  be  found  from  the  56  H*Wind 
plots  (see  Fig.  11).  Generally,  strong  hurricane  conditions  agreed 
with  this  theoretical  expectation  the  of  storm  direction  and  radius 
to  A  being  collinear,  however,  for  other  hurricanes  the  storm 
speed,  storm  direction,  and  central  pressure  affect  the  orientation 
of  A,  B,  and  C.  Using  the  56  plots  we  attempted  to  determine  the 
angles  between  the  radius  to  A  and  the  storm  track  (0C).  Unfortu¬ 
nately,  no  simple  pattern  emerged.  Therefore,  the  data  points  were 
placed  into  6  categories  based  on  storm  direction.  While  no  mean¬ 
ingful  difference  between  the  two  vectors  was  observed  for 
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0  <  0  <  180,  for  the  other  sectors,  the  difference  varied  with  Pc 
and  Vs.  The  results  are  summarized  in  Table  3.  The  net  0  used  in 
Eq.  (7)  is,  therefore,  the  sum  of  0  r  and  0  c. 

Figure  12  shows  plots  for  the  wind  fields  pertaining  to  H*Wind, 
the  RV  model,  the  RV  model  with  the  d(Q)  correction  before 
applying  the  0C  correction,  and  the  RV  model  with  the  corrected 
0.  The  hurricane  shape  in  Fig.  12(d)  is  a  clearly  substantial 
improvement  over  the  basic  RV  model.  However,  the  hurricane 
size  and  the  velocity  distribution  are  not  well  matched;  thus,  fur¬ 
ther  adjustment  is  required. 

3.2  Development  of  Correction  Factor  CB .  Imperfections  at 
any  grid  point  in  the  modeled  velocity  Vrv  1  are  due  to  two  reasons. 
The  first  is  due  to  the  fact  that  the  modeled  velocity  is  not  at  the  10 
m  elevation,  whereas  H*Wind  is;  the  second  is  due  to  errors  in  the 
hurricane  size  and  velocity  distributions.  The  first  error  is  addressed 
in  this  paper  as  a  function  of  r,  and  the  second  as  a  function  of  both 
r  and  0.  In  other  words,  we  provide  a  correction  to  the  Vrv  1  using  a 
correction  factor  CB  that  includes  both  sources  of  errors,  defined  as 

CB{r,  9)  =/i(r)  +/2(r,  9)  (8) 

The  function  /j  (r)  was  chosen  as  a  constant  (  =  0.8)  by  Phadke 
et  al.  [25].  Figure  13  shows  an  example  of  the  RV  modeled  wind 
and  the  corresponding  data  along  the  central  transect  for  hurricane 
Katrina  (0600  UTC,  8/27/2005),  with  an  adjustment  by  a  factor  of 
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Ivan 


Katrina 


Rita 


Dolly 


Fig.  5  Modeled  and  measured  wind  velocities  along  the  central  transect 


0.8  (Fig.  13 (/?)).  It  is  obvious  that  this  adjustment  leads  to  an 
underestimation  outside  the  core.  An  examination  of  several  such 
plots  suggested  that  instead  of  using  a  constant,  the  following 
function  was  to  be  preferred 

/lW=08  +  0.2(^), 

/iW  =  l,  for  r  >  4Rm 
fi{r)  =  0.8,  for  r  <  Rm 

041102-8  /  Vol.  134,  NOVEMBER  2012 


This  leads  to  an  improvement,  as  shown  in  Fig.  13(c). 

As  to  the  correction  necessitated  by  varying  storm  sizes,  the 
function f2(r,  6)  was,  for  simplicity,  described  as 

f2(r,0)=F2(r)-E(0)  (10) 

In  order  to  estimate  F2(r)  and  E(0),  the  discrepancies  A(r,  6) 
between  Vrv i  and  H*Wind  at  approximately  100  points  along 
four  radii  (corresponding  to  0  =  0°,  90°,  180°,  and  270°)  were 
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Day 


ranges.  However,  the  similarity  of  the  curves  in  the  r-direction  for 
different  directions  justifies  the  separation  of  variables  chosen  in 
Eq.  (10).  For  the  storms  with  36  m/s  <  Vm  <  50  m/s,  the  maximum 
discrepancy  for  each  0,  denoted  by  A m(0),  occurred  at  approxi¬ 
mately  r  =  4 Rm.  Figure  14  (right)  shows  a  composite  of  all  dis¬ 
crepancies  for  such  cases.  A  regression  analysis  was  used  to  fit  the 
following  curve  to  rectify  these  errors 


AM) 

AM) 


=  0 


for  r  <  Rm 


Fig.  6  SWH  comparison  plot  during  Hurricane  Katrina  (Aug. 
2005)  at  NDBC  buoy  42  040,  located  approximately  174  km  from 
the  storm  track 


AM) 

AM) 


— 0.0218r3  +0. 1017r2 +0.281/- - 0.3609,  Rm<r< 4 Rm 

(11) 


determined.  From  the  400  values  of  A(r,  6 ),  the  maximum  dis¬ 
crepancy  Amax  was  determined. 

In  general,  two  patterns  for  the  discrepancies  were  observed,  as 
shown  in  Figs.  14  and  15  (left  frames),  based  on  wind-speed 


=  -  0.000333r*  +  0.01639 r2  +  0.12041r  -  0.7565, 

A„(0) 

x  4 Rm  <  r  <  15 Rm 


Fig.  7  H*Wind  (left)  and  RV-model  (right)  wind  speeds  (m/s):  Katrina  (0600  UTC  28  August.  2005),  Pc  =  930  (top);  Ike  (1200  UTC 
10  September.  2008),  Pc  =  959  (bottom) 
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9M  12D  18:00  Wind  Speed(m/s) 


Fig.  8  Deformation  schematic 


Table  2  Deformation  ratios 


Number  of 
HWIND  wind 

/  \ 

field  data  dA 

=  (vA/vB)  dB  = 

{Vb/Vb) 

dc  =  ( Vc/Vb ) 

Ivan 

12 

1.114 

1 

0.842 

Katrina 

11 

1.142 

1 

0.875 

Rita 

13 

1.166 

1 

0.872 

Dolly 

7 

1.180 

1 

0.703 

Ike 

13 

1.160 

1 

0.838 

Average 

1.152 

1 

0.826 

1.2 
1 

0.8 
tj°  0.6 
0.4 
0.2 
0 

925  930  935  940  945  950  955  960  965  970 

Central  Pressure  (mb) 

Fig.  9  Best-fit  line  for  dc  vs  Pc 


For  Vm  >  70  m/s,  the  pattern  is  largely  reversed  and  the  data 
(Fig.  15)  suggested  using  the  following  form 


For  50  m/s  <  Vm  <  70  m/s,  the  discrepancy  between  the  modeled 
wind  speed  and  H*Wind  was  negligible  for  the  most  part. 

Note  that  the  left  hand  side  of  Eqs.  (11)  and  (12)  is  a  measure 
of  F2  (r)  because 

AM)  F2(r)-E(0)  F2(r) 

A m(0)  F2{4Rm)  ■  E(6)  F2(4Rm) 


As  a  result  of  Eqs.  (8)  and  (13),  it  may  be  seen  that  E(0)  would 
represent  the  normalized  velocity  discrepancies  along  a  circle 
with  r  =  4Rm.  To  estimate  E(0),  the  discrepancy  along  this  circle 
at  0  =  0°,  90°,  180°,  and  270°  were  used.  These  four  locations  are 
denoted  by  A,  B,C,  and  D,  and  the  discrepancies  at  these  locations 
are  denoted  by  EA,  EB,  Ec,  and  E°.  A  curve  was  fitted  to  the  val¬ 
ues  at  these  four  locations  to  estimate  E(0)  for  other  values  of  6 , 
in  a  manner  analogous  to  that  for  d(0).  An  examination  of  the  data 
showed  that  these  normalized  discrepancies  at  A,  B,  C,  and  D  varied 
with  Vm.  As  was  the  case  for  the  radial  error,  the  parameters  showed 
different  characteristics  for  large  wind  speeds  (Vm  >  60  m/s)  and 
for  moderate  wind  speeds  (36  m/s  <  Vm  <  50  m/s).  A  detailed  ex¬ 
amination  of  the  discrepancies  led  to  the  following  expressions 


Eb  =  0.4(Vm  -  70)/ 10, 
Eb  =  0.9(50  -  Vm)/20 
Eb  =  0,  otherwise 
Ed  =  0.5(Vm  -  60)/20, 
Ed  =  0.7(50  -  Vm)/20, 
Ed  =  0,  otherwise 


for  Vm  >  70 
for  36  <  Vm  <  50 

for  Vm  >  60 
for  36  <  Vm  <  50 


(14a) 


(14  b) 


^^=4.868e-V  — 0.001871/-3+0.02643>-2— 0.02104r+0.115 

Am(9) 


(12) 


In  general,  we  found  that  EA  «  Ec  «  ( EB  +  E°)/2. 

With  these  parameterizations,  the  function  E(0)  was  defined 
as  E(0)  =  a\  cos(2 6)  +b\  for  0  <  9  <  n  and  E(0)  =  a2  cos(2 9) 


Fig.  1 0  The  d(0)  for  different  central  pressure 
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Fig.  11  Wind  speeds  (m/s)  during  Hurricane  Ike  (2008,  0d=  305°):  H*Wind  (left),  and  RV-model  (right) 


Table  3  Storm  direction  dependence  on  central  pressure  and 
translation  speed 


Sector 

Pc  (mb) 

Vs  (m/s) 

ec 

180  <  6  <  255 

-90 

255  <  6  <  285 

<  3.88 

90 

V,  >  3.88 

90  —  60  (1  — y) 

285  <  6  <  315 

Pc  >  970 

V,  >  3.88 

60  ■  (1  -  y) 

Pc  <  970 

V,  <  3.88 

ON 

O 

>  3.88 

60  ■  a  +  60  •  (1  —  y) 

315  <  6  <  340 

Pc  >  970 

V,  >  3.88 

60  ■  (1  —  y) 

Pc  <  970 

V,  <  3.88 

60  •  a 

Vs  >  3.88 

60  ■  (1  -  y) 

340  <  6  <  360 

>  3.88 

60  •  (1  —  y)/2 

Note:  y  =  (5.55  -  Vs)/1.67  (for  3.88  <  <  5.55),  7  =  0  (for  Vs  >  5.55), 

and  7  =  1  (for  V,  <  3.88). 
aa  is  defined  in  Eq.  (6). 


+  b2  for  7i  <  6  <  271,  where  the  constants  a\,  <22,  b\,  and  £2  are 
determined  from  the  conditions  E(0)  =  ED  at  6  =  3/2n  and 
E(0)  =  EB  at  6  =  7t/2.  The  final  form  of  E(0)  reduced  to 

E(0)  =  [(. Ed  -  Eb)  •  cos  2  6  +  (3EB  +£d)]/40  <6  <n 
E{9)  =  \{Eb  -Ed)  ■  cos  29  +  ( Eb  +  3E°)}/4n<  9  <  2n  (15) 

To  conclude,  a  detailed  examination  of  five  storms  (approximately 
56  plots)  suggested  that  the  10  m  elevation  wind  could  be 
obtained  by  adjusting  the  RV  model  as  follows 

Vcrv  =  Vrv  •  d(0)  •  Cb  (16) 

The  effect  of  this  second  adjustment  may  be  seen  in  Figs.  12(e) 
and  12(f)  for  Hurricane  Ike  and  in  Fig.  16  for  Hurricane  Katrina. 
We  note  that  in  the  case  of  Hurricane  Katrina,  the  shape  obtained 
by  the  original  RV-model  is  reasonable,  although  for  both  storms, 
the  final  adjusted  velocities  based  on  eq.  (15)  are  significantly 
improved  as  compared  to  the  original  RV  model  results.  The 
effect  of  our  modifications  on  the  wave  heights  during  Hurricane 
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Katrina  is  shown  in  Fig.  6;  the  differences  are  of  a  magnitude  suf¬ 
ficient  to  influence  extreme  wave  height  calculations. 

4  Validation  and  Effect  on  Wave  Fields 

As  noted  earlier,  the  Reanalysis  wind  fields  are  available  on  a 
6-hourly  basis,  dating  back  to  1948,  and  have  occasionally  been 
used  for  wave  and  storm-surge  hindcasting  [16-20,37].  However, 
the  coarse  resolution  (2.5°  x  2.5°)  can  be  expected  to  adversely 
impact  the  wind  and  wave  hindcasts  in  hurricane  conditions.  For 
example,  considering  the  case  of  Hurricane  Gordon,  as  mentioned 
in  Sec.  1,  winds  resulting  from  the  modified  RV  model  are  shown 
in  Fig.  17.  The  wave  hindcast  resulting  from  these  wind  fields  is 
compared  with  data  from  buoy  42003  in  Fig.  18.  It  is  clear  that 
the  Reanalysis  winds  are  inadequate  for  reliable  wave  height  pre¬ 
dictions  and  that  the  modified  RV  model  results  in  better  wave 
predictions,  with  a  maximum  difference  on  the  order  of  4  m. 

Nonetheless,  the  Reanalysis  winds  can  be  used  to  provide  the 
“background”  windfields,  viz.,  the  windfields  outside  the  hurri¬ 
cane  and  also  the  windfields  that  existed  before  the  onset  of  the 
hurricane.  The  benefits  of  incorporating  such  background  wind- 
fields  have  been  indicated  by  Liu  et  al.  [35].  We  have,  therefore, 
constructed  a  composite  windfield  dataset  by  merging  the  modi¬ 
fied  RV  winds  with  the  Reanalysis  winds  (the  latter  are  interpo¬ 
lated  onto  a  5  km  grid). 

The  parametric  model  equation  (Eq.  (16))  was  tested  against  17 
storms  (107  plots)  covering  a  period  from  1961  to  2008  (Table  4); 
these  formed  a  subset  of  the  entire  H*Wind  dataset.  The  others 
were  not  chosen  because  numerical  data  were  not  available  for 
them  or  because  they  contained  too  many  land  points. 

For  each  snapshot  V(x,y)  was  calculated  using  Eq.  (15)  for  a 
matrix  containing  81  x  81  points  on  each  side  (at  a  resolution  of 
approximately  5  km)  covering  a  4°  x  4°  region  (approximately 
6561  grid  points  per  snapshot).  Wind  speeds  from  the  modified 
RV  model  and  the  original  RV  models  were  compared  with  the 
H*Wind  data  for  these  grid  points  for  each  snapshot.  The  average 
of  all  of  the  absolute  values  of  the  percentage  error  at  each  grid 
point,  denoted  by  AE,  is  presented  in  Table  5.  In  general,  in  14 
out  of  17  cases,  the  modified  RV  model  (Vcrv)  produces  far  lower 
errors  than  the  original  model.  Disregarding  small  differences 
(say,  less  than  10%,  which  occur  in  eight  cases),  eight  of  the 
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Fig.  12  Wind  speeds  (m/s)  for  Hurricane  Ike  at  UTC  1200  11  Sept.  2008;  0d=  295°,  Vs=  4.44  m/s,  Pc=  946  mb  (a)  H*Wind,  ( b ) 
RV-model,  (c)  before  angle  correction  with  d(0),  (d)  after  angle  correction  with  d(0),  (e)  after  application  of  both  correction  fac¬ 
tors,  and  (/)  comparison  along  the  central  E-W  transect 
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Fig.  13  Wind  speed  along  the  central  transect  for  the  (a)  H*Wind  and  RV  model,  ( b )  H*Wind  and  RV-model  with  10-meter  eleva¬ 
tion  correction  factor  (  =  0.8),  and  (c)  H*Wind  and  RV-model  with  the  10-meter  elevation  correction  factor  (from  Eq.  (9)) 


Fig.  14  Model-data  discrepancies  for  Hurricane  Ike:  Vm  =  48.6  m/s  (left),  and  best-fit  curve  based  on  400  data  points  (right) 


Fig.  15  Model-data  discrepancies  for  Hurricane  Rita  with  Vm  =  79.2  m/s  (left),  and  best-fit  curve  based  on  400  data  points 
(right) 


remaining  nine  showed  significant  improvements,  while  only  one 
(i.e.,  Hurricane  Dennis)  showed  considerable  deterioration  rela¬ 
tive  to  the  original  RV  model.  The  maximum  improvement  is  for 
Hurricane  Dean  (30%).  Unfortunately,  for  the  two  cases  with 
maximum  discrepancies  (i.e.,  Hurricanes  Dennis  and  Dean),  the 
assessment  depends  on  very  few  snapshots  being  available.  In 
addition,  a  detailed  comparison  of  both  parametric  wind  fields 
against  the  H*Wind  data  was  conducted  for  Hurricane  Georges 
(1998)  and  the  results  are  in  shown  in  Fig.  19. 

Journal  of  Offshore  Mechanics  and  Arctic  Engineering 


The  results  indicate  clearly  that  the  CRV-model  wind  speed 
contour  shape  shows  a  better  match  with  the  H*Wind  than  the 
RV-model  plot.  The  wave  fields  resulting  from  the  two  numeri¬ 
cally  generated  wind  fields  are  shown  in  Fig.  19  (bottom).  For 
both  wind  and  wave  fields,  although  the  maximum  values  are 
largely  the  same,  the  spatial  patterns  are  different,  leading  to  a 
maximum  of  about  a  ±1.8  m  difference  in  wave  height  estimates 
at  some  locations.  A  difference  on  the  order  of  about  3.5  m  in  the 
maximum  SWH  is  obtained  for  Hurricane  Katrina  (see  Fig.  6 
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Fig.  16  Wind  speeds  (m/s)  during  Hurricane  Katrina  (0600  UTC  28  Aug.  2005);  (a)  H*Wind,  ( b )  RV-model,  (c)  CRV-model,  and 
( d)  comparison  along  the  central  transect 


Fig.  17  Wind  speed  during  Sept.  2000  at  NDBC  buoy  42  003 


where  buoy  data  are  also  shown)  and  about  4  m  for  Hurricane  Ike 
(not  shown).  Such  large  differences  have  the  potential  to  affect 
the  calculations  of  extreme  wave  statistics. 

Finally,  the  effect  of  the  two  wind  fields  ( VRV  and  VCrv )  on  the 
100-yr  return  period  SWHs  (defined  by  SWH100)  are  examined. 
These  extreme  wave  height  calculations  were  made  by  fitting  the 
Gumbel  distribution  to  51  years  of  simulations  (1958-2008).  The 
purpose  of  this  exercise  was  merely  to  study  the  differences; 
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Fig.  18  SWHs  during  Hurricane  Gordon  (Sept.  2000)  at  NDBC 
buoy  42003  (located  approximately  57  km  from  storm  track) 


actual  estimates  of  SWH100  would  depend  on  the  type  of  statisti¬ 
cal  model  used,  and  other  factors  which  will  be  discussed  in  a  sep¬ 
arate  paper.  The  results  (Fig.  20)  show  considerable  differences  in 
the  spatial  patterns.  These  differences  (on  the  order  of  4  m  at 
some  locations;  see  Fig.  20(c)),  are  substantial  relative  to  the 
SWH10o  obtained  using  the  original  RV  model  and  also  relative  to 
the  values  used  by  API  (2000)  [1]. 
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Table  4  Hurricanes  selected  to  test  Eq.  (1 6) 


Table  5  Performance  of  RV  and  CRV  models 


Number  of 

AE  (CRV  model) 

AE  (RV  model) 

HWIND  storm 

No. 

Hurricane 

N 

(%) 

(%) 

No. 

Name 

Storm  period 

snapshots  chosen 

1 

Carla 

2 

15 

20 

1 

Carla  (1961) 

Sep.  03.  12:00-Sep.  16.  00:00 

2 

2 

Georges 

9 

12 

26 

2 

Georges  (1998) 

Sep.  15.  12:00-Oct.  01.  06:00 

9 

3 

Gordon 

5 

22 

26 

3 

Gordon  (2000) 

Sep.  14.  12:00-Sep.  21.  06:00 

5 

4 

Isidore 

5 

37 

45 

4 

Isidore  (2002) 

Sep.  14.  18:00-Sep.  27.  18:00 

5 

5 

Lili 

4 

23 

16 

5 

Lili  (2002) 

Sep.  21.  18:00-Oct.  04.  12:00 

4 

6 

Ivan 

12 

11 

24 

6 

Ivan  (2004) 

Sep.  02.  18:00-Sep.  24.  06:00 

12 

7 

Arlene 

4 

32 

34 

7 

Arlene  (2005) 

June  08.  18:00-June  14.  06:00 

4 

8 

Dennis 

3 

38 

20 

8 

Dennis  (2005) 

July.  04.  18:00-July  18.  06:00 

3 

9 

Emily 

6 

26 

20 

9 

Emily  (2005) 

July  11.  00:00-July  21.  12:00 

6 

10 

Katrina 

9 

9 

20 

10 

Katrina  (2005) 

Aug.  23.  18:00- Aug.  31.  06:00 

9 

11 

Rita 

10 

10 

16 

11 

Rita  (2005) 

Sep.  18.  00:00-Sep.  26.  06:00 

10 

12 

Wilma 

12 

11 

30 

12 

Wilma  (2005) 

Oct.  15.  18:00-Oct.  26.  18:00 

12 

13 

Alberto 

3 

19 

24 

13 

Alberto  (2006) 

June  10.  06:00-June  19.  06:00 

3 

14 

Dean 

2 

12 

42 

14 

Dean  (2007) 

Aug.  13.  06: 00- Aug.  23.  00:00 

2 

15 

Dolly 

4 

20 

33 

15 

Dolly  (2008) 

July  20.  12:00-July  27.  00:00 

4 

16 

Gustav 

5 

8 

22 

16 

Gustav  (2008) 

Aug.  25.  00:00-Sep.  05.  12:00 

5 

17 

Ike 

12 

14 

39 

17 

Ike  (2008) 

Sep.  01.  06:00-Sep.  15.  12:00 

12 

Longitude  (W)  Longitude  (W)  Longitude  (W) 


Fig.  19  Wind  speeds  (m/s)  and  SWHs  during  Hurricane  Georges  (at  0600  UTC  27  Sept.  1998):  (a)  H*Wind,  ( b )  CRV-model  wind, 
(c)  RV-model  wind,  ( d)  SWH  (m)  using  CRV-model  wind,  (e)  SWH  (m)  using  RV-model  wind,  and  ( /)  the  difference  SWH  (m) 
between  ( d)  and  (e) 


5  Conclusions 

An  integrated  set  of  offshore  and  coastal  models  for  winds, 
tides,  waves,  and  storms  (hurricanes  and  typhoons)  is  extremely 
useful  in  the  planning  and  in  the  management  of  an  emergency 
response  pertaining  to  tropical  storms  that  affect  US  coasts  and 
islands.  These  models  are  also  used  in  other  areas  of  coastal  engi¬ 
neering  concerned  with  the  design  of  coastal  and  wetland  protec- 
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tion  systems  such  as  barrier  islands,  levees,  and  coastal  navigation 
and  harbor  structures.  In  these  applications,  users  need  models 
with  varying  complexity  and  accuracy  for  calculating  winds  when 
modeling  past  storm  events  (hindcast),  future  storms  (forecast), 
and  hypothetical  (synthetic)  storms.  Less  resource-demanding 
models  that  run  fast  with  a  low-order  accuracy  may  be  used  to 
perform  quick  simulations  for  decision-making  and  to  determine 
if  further  modeling  using  models  of  higher  accuracy  is  warranted. 
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Fig.  20  SWH100  (m)  for  the  Gulf  of  Mexico  estimated  using  (a)  RV-model  wind,  ( b )  CRV-model  wind,  and  (c)  differences  between 
(a)  and  ( b ) 


The  parametric  wind  modeling  approach  described  in  this  paper  is 
intended  to  address  these  critical  needs,  which  can  be  used  to  esti¬ 
mate  winds  when  no  other  information  is  available  for  historical 
storms  of  interest. 

We  have  investigated  the  RV  model  and  other  similar  recent 
parametric  wind  models  with  corrections  for  a  large  number  of 
tropical  storms  in  the  GOM.  We  found  that  these  models  produced 
a  substantial  difference  in  wind  speeds  relative  to  the  H*Wind 
data  in  the  GOM.  This  motivated  us  to  consider  spatial  adjust¬ 
ments  to  the  shape  and  distributions  to  the  contours  produced  by 
the  RV  model.  We  developed  an  adjusted  model  (CRV)  applica¬ 
ble  to  asymmetric  wind  fields  typically  observed  in  tropical 
storms,  and  have  conducted  a  thorough  testing  of  it  against 
H*Wind  for  several  hurricanes. 

The  CRV  model,  for  the  most  part,  produced  a  better  shape  and 
distribution  of  the  velocity  plots.  In  a  comparison  performed  for 
17  hurricanes  with  56  snapshots,  the  CRV  model  produced 
smaller  percentage  errors  than  the  RV  model.  The  CRV  wind 
fields  produced  SWHs  which  were  on  the  order  of  4  m  higher  for 
some  hurricanes  examined.  Results  presented  in  Figs.  6  and  18 
show  a  better  match  with  data.  The  estimates  of  SWH100  resulted 
in  a  4-m  difference  (maximum)  at  some  locations. 
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